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ABSTRACT 

We have reanalyzed the 260 ksec XMM-Newton observation of 1E1207. 4-5209. There are several 
significant improvements over previous work. Firstly, a much broader range of physically plausible 
spectral models was used. Secondly, we have utilized a more rigorous statistical analysis. The standard 
F-distribution was not employed, but rather the exact finite statistics F-distribution was determined by 
Monte-Carlo simulations. This approach was motivated by the recent work of Protassov et al. (2002) and 
Freeman et al. (1999). They demonstrated that the standard F-distribution is not even asymptotically 
correct when applied to assess the signficance of additional absorption features in a spectrum. With our 
improved analysis we do not find a third and fourth spectral feature (Bignami et al. (2003), De Luca et 
al. (2004)) in 1E1207. 4-5209, but only the two broad absorption features previously reported (Sanwal 
et al. (2002), Hailey and Mori (2002)). Two additional statistical tests, one line model dependent and 
the other line model independent confirmed our modified F-test analysis. For all physically plausible 
continuum models where the weak residuals are strong enough to fit, the residuals occur at the instrument 
Au-M edge. As a sanity check we confirmed that the residuals are consistent in strength and position 
with the instrument Au-M residuals observed in 3C273. 

Subject headings: stars: neutron - X-rays - individual: 1E1207. 4-5209 - methods: statistical 



1. INTRODUCTION 

The thermal emission from isolated neutron stars (INS) 
probes neutron star physics since it is uncontaminated by 
emission from an accretion disk. In particular, many INS 
(age ~ 10 3 -10 5 yrs) are good candidates because they are 
hot enough that the thermal emission is not obscured by 
non-thermal magnetospheric emission. Lines in thermal 
spectra are an excellent tool to probe the strong magnetic 
field and gravity of the INS. 

Two broad absorption features were discovered in the 
INS 1E1207 by Chandra (Sanwal et al. 2002) and later 
confirmed by XMM-Newton (Mereghetti et al. 2002). The 
discovery was followed by the detection of a single spectral 
feature in the X-ray thermal spectra of three nearby INS 
(Haberl et al. 2003; van Kerkwijk et al. 2004; Haberl et al. 
2004). 1E1207 still remains unique because it shows more 
than one spectral feature as opposed to the other three 
INS which only show a single absorption feature. 

The interpretation of features in NS thermal spectra is 
not straightforward because of effects such as magnetic 
field, gravity and unknown surface composition. Soon af- 
ter the discovery of features in the 1E1207 spectrum sev- 
eral different interpretations were proposed ranging from 
Helium atomic lines at B ~ 10 14 G (Sanwal et al. 2002), 
Iron atomic lines at B ~ 10 12 G (Mereghetti et al. 2002), 
Oxygen/Neon atomic lines at B ~ 10 12 G (Hailey & Mori 
2002; Mori & Hailey 2003) and hydrogen molecule ion lines 
at B ~ 10 14 G (Turbiner & Lopez Vieyra 2004). Recently, 
Bignami et al. (2003) (hereafter B03) proposed electron 
cyclotron lines based on the detection of two additional 
absorption features at 2.1 and 2.8 keV in the 260 ksec 
XMM-Newton observations. They interpreted these four 



absorption lines as the fundamental and three harmonics 
of cyclotron lines. De Luca et al. (2004) (hereafter DL04) 
re-analyzed the same XMM-Newton data and also found 
that the 3rd and 4th feature are significant. 

In this paper we present our detailed analysis of the 
same 260 ksec XMM-Newton observation and investigate 
the significance of the two additional absorption features 
at 2.1 keV and 2.8 keV. 

2. DATA REDUCTION METHODOLOGY 

XMM-Newton observed the INS 1E1207 for two orbits 
starting on August 4, 2002. The PN camera was operated 
with a thin filter in the "small window" mode. This allows 
accurate timing (6 ms time resolution) and minimizes pile 
up. The two MOS cameras were operated in "full frame" 
mode with a thin filter. Our data reduction followed the 
latest XMM-Newton calibration document (Kirsch 2004) 
and we processed the data by running the pipeline directly 
from the ODF files with SAS version 5.4.1. 3 We used 
canned response matrices for all the instruments and gen- 
erated ancillary files by the SAS command " arfgen" . We 
also generated response function using the SAS command 
"rmfgen" , but the difference in chi-squared and fit param- 
eters was tiny. Following Kirsch (2004), we generated a 
lightcurve of high energy (> 10 keV) single pixel (PAT- 
TERN = 0) events on the full field of view and selected 
good time intervals using the criteria for count rates < 1.0 
cts/sec (PN) and 0.35 cts/sec (MOS). The total exposure 
time for PN and MOS after removing time intervals with 
high background are 156 ksec and 231 ksec respectively. 

We extracted the source photons from a 45" radius circle 
centered on the source for all cameras. There arc sufficient 
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photons in the PN camera to enable selection of only sin- 
gle events to achieve the best spectral resolution. Most of 
our results are based on this higher quality singles data, 
however later we discuss some results from the doubles 
analysis. The MOS cameras have less photons so we used 
all events in our analysis to improve the photon statistics. 

Background regions were selected according to the 
recommendations of the XMM-Newton calibration team 
(Kirsch 2004). In the PN camera we chose a 45" ra- 
dius circle background region the same distance from 
the readout node as the source region to ensure simi- 
lar noise. An annular background region encircling the 
source was not used to avoid out-of-time events from the 
source. In the MOS cameras we used a circular region 
away from the source according to the calibration team 
recommendations. We also tried an annular region sur- 
rounding the source and found no significant deviations 
in the continuum fit parameters. The source is located 
in a supernova remnant and non-uniform X-ray emis- 
sion from the remnant could influence the neutron star 
spectra. To validate our choice of background region we 
chose three different background regions for each instru- 
ment and confirmed the spectral fit parameters of the 
source did not significantly change with our choices. As 
a result, source count rates in 0.2-5.0 keV are 1.546+/- 
0.003 cts/sec (PN single+double), 0.346+/-0.001 cts/sec 
(MOS1) and 0.353+/-0.001 cts/sec (MOS2). Background 
count rates in the energy band are 0.089+/-0.001 cts/sec 
(PN single+double), 0.0065+/-0.0005 cts/sec (MOS1) and 
0.0069+/-0.0005 cts/sec (MOS2). 

The spectra were rebinned such that each bin contained 
at least 40 counts and analyzed using XSPEC vl 1.3.1 in 
the 0.3-4.0 keV energy range. Modest oversampling of the 
energy resolution kernel was utilized. Our number of bins 
and degrees of freedom are comparable to those of other 
analyses of this source. We adopted the energy band be- 
tween 0.3 and 4.0 keV for the PN instrument. For the MOS 
instruments we selected the 0.5-4.0 keV band because the 
MOS is not well-calibrated at lower energies. Continuum 
fit parameters for the PN and MOS instruments were in 
good agreement. 

The long XMM-Newton observation permits splitting 
the full observation in order to search for time variabil- 
ity. We split the observation into four approximately equal 
time intervals and found no deviations between the four 
observations. They yielded statistically consistent contin- 
uum fit parameters. We conclude that there was no time 
dependence during the total observation and therefore we 
use the data from the entire observation for all subsequent 
analysis. 

3. ANALYSIS OF DATA USING CONTINUUM MODELS 

3.1. Introduction 

In this section we consider fits to the global 1E1207 spec- 
tra using a variety of physically motivated spectral mod- 
els. Most of these models are not compatible with the data 
and are rejected. For those models which provide good fits 
to the data we determine the relevant fitting parameters. 
The issue of properly fitting additional spectral features 
is a complex one, and we postpone addressing it until §4. 
Here we only set the stage by providing information on the 
best fit continuum models. 



3.2. Atmosphere models for neutron star thermal spectra 

The fitted continuum models must reflect the full range 
of plausible INS atmospheric conditions. The surface of an 
INS can be any element from Hydrogen to Iron, since only 
~ 1O~ 19 M of material produces an optically thick atmo- 
sphere (Romani 1987). The atmosphere distorts the spec- 
tral energy distribution (SED) from Planckian, as will the 
presence of a magnetic field. For B = 0, hydrogen atmo- 
sphere models have the hardest spectrum due to free-free 
absorption. INS spectra soften and show a proliferation 
of lines and edges at higher energies, becoming closer to 
black-body (Romani 1987). For increasing B and atomic 
number there is significant spectral softening of the SED 
(Ho & Lai 2003). Even for Hydrogen, electron binding 
effects can amount to an appreciable fraction of the INS 
surface kT at high enough B-field, leading to incomplete 
ionization and softening of the SED due to opacity from 
bound species (Ho et al. 2003). From the above consider- 
ations we conclude that realistic atmospheric models for 
1E1207 will have an SED somewhere between black-body 
and a fully-ionized hydrogen model, and our subsequent 
fitting reflects this range of SED. 

3.3. Description of model fits to the data (no third or 
fourth line assumed) 

In the following few sections we describe our fit to the 
data without the assumption of a third or fourth line in 
the spectrum. Fits with such additional lines are discussed 
in §4. Therefore our "baseline" models always consist of 
the two universally accepted lines in the 1E1207 spectrum, 
which we model as Gaussian absorption lines at 0.7 and 
1.4 keV. 

For each continuum component we adopt one of the 
following three physically plausible models; blackbody 
(BB), magnetized hydrogen atmosphere (HA) or power- 
law (PL). The magnetized hydrogen atmosphere model 
was constructed assuming a fully-ionized hydrogen atmo- 
sphere at B = 10 12 G, close to the dipole field strength 
as measured from the spin-down parameters (private com- 
munication with V.E. Zavlin). We also tried B-fields as 
low as B = (Zavlin et al. 1996), and this did not alter 
the continuum parameters significantly. Higher B-ficlds 
( ^4.4 x 10 13 G), just give results closer to the black-body 
case (Ho & Lai 2003). 

In our first fits to the spectrum we used one continuum 
component along with the two absorption lines. For all 
instruments the one-component models did not yield ac- 
ceptable x 2 values (xl > 1-3) and showed significant resid- 
uals above ~ 2 keV. This suggested the need for another 
continuum component. 

Given two continuum components and two absorption 
lines there are three types of continuum models defined as 
below. They consist of two continuum components (CI 
and C2) and two absorption lines (LI and L2). CI /LI 
and C2/L2 refer to the lower and higher energy spectral 
component respectively. The three classes of models con- 
sidered (with the notation C * L meaning line L resides on 
continuum component C) 

Model I: CI * LI * L2 + C2 

Model II: CI * LI + C2 * L2 

Model III: (CI + C2) * LI * L2 
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For LI and L2 the Gaussian absorption line was modeled 

as 

F(E) oc exp{~Texp(-(S - E ) 2 /2w 2 )} (1) 

where r and w refer to the line depth and width. We 
fit other absorption line profiles such as a Lorentzian but 
the results were similar. While the exact shape of the ab- 
sorption features may be asymmetric or have substructure 
from blended lines (Mori & Hailey 2003), we obtained ex- 
cellent fits without invoking asymmetric profiles. Thus we 
chose simple Gaussian line shapes. When we fit photo- 
absorption edges to the two absorption features at 0.7 and 
1.4 keV, the xt was n °t as good as for the Gaussian lines 
in any of the continuum models. 

3.3.1. Physical meaning of model I, II and III 

Interpretation of spectral features in NS thermal spectra 
is not straightforward since the line parameters depend on 
various NS parameters such as surface element, ionization 
state, magnetic field strength and gravitational redshift 
(Hailey & Mori 2002; Mori & Hailey 2003). Therefore, it 
is important whether the observed spectral features origi- 
nate from the same region or layer. 

The presence of two continuum components complicates 
our interpretation of the observed absorption features. We 
assume that the two continuum components originate from 
different regions on the NS surface. For instance, CI is 
emitted from a large area on the surface and C2 is emit- 
ted from a hot polar cap (DL04). 

Model I assumes that LI and L2 are from the same re- 
gion as the emission of CI and Model II assumes that they 
are from different regions. Model III assumes that LI and 
L2 are from a layer above the two regions emitting contin- 
uum photons. Several physical models predicting a layer 
made of electron-positron or electron-ion pairs a few NS 
radii above the NS surface have been proposed (Dermer & 
Sturner 1991; Wang et al. 1998; Ruderman 2003). Such a 
layer may become optically thick and modify the thermal 
spectra from the NS surface. In addition, absorption in 
model III can take place anywhere between the NS and the 
observer (e.g. magnetosphere, supernova remnant, ISM or 
materials on the XMM-Newton mirrors and cameras). 

3.4. Results of continuum model fitting to XMM-Newton 

data 

We searched for continuum models that are consistent 
with the data using models I, II and III. The set of models 
we considered is extensive including all permutations of 
blackbody (BB), magnetized hydrogen atmosphere (HA) 
and power law (PL) models for CI and C2. 

All the two-component thermal models (and combina- 
tions of BB and HA) fit the data well, yielding xl ver Y 
close to unity (table 1 and 2). For brevity we do not 
consider mixed cases (BB+HA) because tests showed that 
in assessing absorption line significances (§4) they always 
produced results intermediate between the pure BB and 
pure HA cases. Continuum models with PL components 
are ruled out because they do not adequately fit the data, 
leaving significant residuals above 2 keV. Models I, II and 
III with thermal continuum components are all acceptable. 

Figure 1 show the spectra and residuals for the best fit 
to model I(BB+BB) using the PN single/double events 



and MOS1/MOS2 data. We also fit the data with a third 
spectral line, calculated x 2 , and evaluated the difference 
of x 2 between third line case and no third line case. This 
is a key parameter in our subsequent statistical tests. 

We show our best-fit continuum parameters and 90% 
confidence levels for six different models in the energy 
range from 0.3-4.0 keV for PN in Table 1. All the er- 
rors quoted are statistical errors and we did not include 
any systematic error. We calculated 90% confidence level 
errors using S = S m i„ + 2.706 for each parameter while al- 
lowing the other parameters to vary freely (Lampton et al. 
1976). These models yield excellent fits to the data with 
X 2 values for PN slightly less than unity. We found that 
Gaussian line profiles yield a statistically good fit to the 
data. The PN blackbody models I, II and III all give con- 
sistent values for these parameters. 

We note that for some of our models the line depths 
(especially for the 1.4 keV line) are not well-constrained. 
This is caused by the presence of two continuum compo- 
nents. If the absorption feature is fixed to reside on one of 
the continuum components then the level of that compo- 
nent becomes very important. If the depth of the absorp- 
tion feature is close to the continuum level the residuals 
cannot be adequately fit rendering the line depth poorly 
constrained. It is possible for both of the broad absorp- 
tion features to lie on the higher temperature continuum 
component (CI + C2 * LI * L2). However, the fit was not 
acceptable leaving significant deviations from the data. 

As an aside we note that we took an additional step to 
confirm the robustness of our continuum parameter fitting 
and x 2 determination. We examined a range of degrees of 
freedom in the PN analysis spanning those previously re- 
ported in the literature, and consistent with obtaining nor- 
mal statistics and sensible energy resolution kernel sam- 
pling. Both our continuum parameter determination and 
quality of fit were insensate to this. However we empha- 
size the extremely important technical point that we se- 
lected the actual bin size a priori, not by adjusting the size 
emipirically to obtain the lowest x 2 ■ It is well known this 
would result in a fit statistic which is not x 2 distributed 
(Eadie et al. 1983). It is important for the subsequent sta- 
tistical analysis that the conditions for a distribution-free 
statistic are met, even in the finite-N limit. 

3.5. Comparison with and improvements to previous 

work 

In several previous investigations (DL04, B03) the data 
were fit with two black-body components and two absorp- 
tion lines. In both cases xl significantly greater than one 
were obtained, requiring the inclusion of additional ab- 
sorption lines to fit some higher energy residuals. Other 
work did not require such extra lines when fitting either 
black-body continuua (Sanwal et al. 2002; Hailey & Mori 
2002; Mereghetti et al. 2002) or magnetized HA models 
(Mereghetti et al. 2002). 

In order to more thoroughly investigate these deviations 
in fit quality using just continuum models we have ex- 
ploited, in addition to more up-to-date response matrices 
and a more time-consuming study of background effects, 
a much broader set of continuum models (models I and II 
in addition to model III). In the case of model III, where 
our comparison with previous work is direct, we find agree- 
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mcnt at the 90% confidence level for all derived parameters 
for both MOS and PN data. 

In regards the overall fits, we find acceptable X 2 (~ 1— 
1.3) for all three classes of models, i.e., our fits do not 
require additional spectral features. However introduction 
of extra lines does improve the overall x 2 very slightly. As- 
sessing the necessity of such additional features is a matter 
of great statistical delicacy, which we discuss fully in §4. 

To facilitate a direct comparison with previous work we 
have used the results of table 2 in DL04 which indicates 
the best continuum model fits (without 2 additional lines) 
as giving xl ~l-8, 1.3 and 1.6 for PN, MOS1 and MOS2 
respectively. While the MOS fits of DL04 are not as good 
as ours (for any of our three models), the discrepancy is 
not really that large; DL04 only marginally required extra 
lines in the MOS data. Only the PN fits appear to have a 
large " discrepancy" . One must be very cautious in assum- 
ing this discrepancy arises from differences in processing of 
the data. Such an assumption is belied by the agreement 
in overall signal and background count rates in MOS and 
PN, in MOS and PN fitting parameters and even in overall 
X 2 for MOS1/2 data between all groups. The problem is 
that in the case of finite count statistics, there can be sub- 
stantial deviations in x 2 from its ideal, distribution-free 
theoretical form (Eadie et al. 1983). 

Simulations support this assertion. We generated 1000 



sets of bootstrapped data based on the 1E1207 data set, 
fitting this data set with continuum model III, and then 
comparing the theoretical and observed x 2 distributions. 
It was found that xl OI 1-6 or greater were found more 
than an order of magnitude more times than predicted 
by the standard x 2 distribution. Thus there is a non- 
negligible probability (<~ 10%) of obtaining a deceivingly 
large x 2 using proper fitting procedures. The difference in 
X 2 for the PN data may be nothing more than artifacts 
of finite N statistics. Our results are more immune to this 
because we fit so many test cases. Moreover such a large 
X 2 is not problematic per se. Rather what is required is 
an appropriate finite-N statistical methodology applicable 
to each reduced data set. An appropriate approach has 
recently been developed and we present it and apply it to 
the 1E1207 data set in §4 and 5. 

4. STATISTICAL TESTS FOR EXTRA SPECTRAL FEATURES 

While it may naively appear that assessing the existence 
of and need for an extra spectral feature in a continuum 
spectrum is straightforward, this is not true. We present 
three approaches for assessing the need for extra spectral 
features. These include goodness-of-fit test for lines in a 
continuum, direct spectral fitting and significance assess- 
ment, and the modified F-test. 

We have chosen to apply the F-test because it has been 
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used in previous work on 1E1207. Despite its great famil- 
iarity to astrophysicists, the use of the F-test to determine 
the need for extra spectral features is extremely problem- 
atic. To effect this analysis we exploit the recent work 
of Protassov et al. (2002) (hereafter P02) and Freeman 
et al. (1999) (hereafter F99). These authors have de- 
veloped approaches which surmount the difficulties in the 
standard F-test and permit much more reliable and accu- 
rate determination of feature significance than is possible 
with the "standard" F-test. Our analysis exactly follows 
that of P02, so we emphasize P02 in our brief discussion of 
theoretical and practical particulars in §4.1. The reader is 
directed to the above references for a more thorough treat- 
ment. In §4.2 and §4.3 we consider the other methods for 
attacking this problem. 

4.1. F-test for detecting extra model components 

In the F-test and its closely related cousin the likelihood 
ratio test (LRT) a test statistic is formed from the distri- 
bution function f(xi\9) governing the probability of counts 
Xi in bin i under a set of parameters 9 — (9\....9 P ). The 
test statistic is a functional containing the f(xi\9) with 9 
free to assume any value and also f(xi\9) where some pa- 
rameters are constrained by the values that are required 
by the null hypothesis. To be concrete, in the case of the 
F-statistic the statistic F — {xi{x, 8) — xK x i #)) Ixt (Bev- 
ington 1969; P02). 

The F-test is extremely powerful because it is 
"distribution-free" . That is, in the asymptotic (large 
count) limit, and under the null hypothesis (specifically, 
in our case, r = 0) it does not depend on the initial prob- 
ability distribution, which is the underlying distribution 
f{xi\9\...9 p ). Rather the test statistics only depend on 
the reference distribution (in our case the F-distribution) . 
The asymptotic behavior effectively decouples details such 
as the shape of the continuum and the particulars of the 
line profile and strength from the way the test statistic 
is distributed asymptotically. This is why the test is so 
powerful and broadly applicable. 

P02 and F99 made the crucial observation that there 
are certain regularity and topological conditions that must 
be met for distribution-free statistics. The regularity 
conditions are related to the integrals and derivatives of 
f{xi\9\...9 p ) (see Serfling 1980; Roe 1992 for a precise 
statement) and are of no interest to us since these con- 
ditions are weak and thus almost always satisfied in astro- 
physics modeling. Of profound importance are the topo- 
logical conditions (P02). These arc: 1) the model under 
test must be nested, that is, the parameter under test in 
the constrained model must be a subset of that parameter 
in the unconstrained model and 2) the parameter under 
test in the constrained model cannot lie on the boundary 
of the possible values of the parameter under the uncon- 
strained test. 

Evaluating these topological conditions for equation 1 
we see that the first is met; the constrained model has 
r = and thus the parameterization in the constrained 
model is a subset of that in the unconstrained model. How- 
ever condition 2 is violated, t = is on the boundary of 
the possible values r can assume under the unconstrained 
model. When condition 2 is violated the F-statistic does 
not even asymptotically approach the F-distribution. P02 



illustrate, through simulations for an absorption line on a 
continuum, that the reference distribution overestimates 
the strength of the line by almost an order of magnitude. 
P02 properly note that in general nothing can be said 
about the significance of a line with the standard F-test, 
because we do not even in principle know the asymptotic 
reference distribution. In this regard, for instance, the 
XSPEC users manual specifically advises against the use 
of the F-test in assessing the necessity for extra lines in a 
spectrum. 

4.1.1. Analysis of 1E1207 line significance using modified 

F-test 

Following P02 we seek to more accurately assess the 
significance of the alleged third and fourth lines in 1E1207 
by employing the F-statistic, which we now recognize is 
not F-distributed. The F-statistic is actually a less obvi- 
ous choice than LRT, but it facilitates direct comparison 
with previous work. Since an appropriate asymptotic dis- 
tribution does not exist, even in principle, what is now 
required is the establishment of a reference distribution 
for finite statistics. P02 and F99 explain in detail how to 
do this. We have determined the reference distribution us- 
ing the posterior predictive p- value methodology (PPPM) 
of P02, and the reader is referred there for details. The 
idea behind PPPM is straightforward although the com- 
puter resources required to implement it are substantial. 
We only outline the procedure here. 

Under the null hypothesis we can use data to generate 
fake data from which we can determine the relevant ref- 
erence distribution. However this standard bootstrapping 
technique is not correct, since we do not know the "true" 
values of the parameters to use in the simulation. The fail- 
ure of the topological condition means we lose conditioning 
properties which would make it meaningful to replace pa- 
rameters by best estimators. The solution, as pointed out 
by P02, is to use parameter values in the simulation which 
are likely given the observed data. This simply requires 
a statistical distribution of parameter values which can 
be sampled and incorporated into the simulation which 
generates the distribution of the statistic. P02 advocate 
a Bayesian approach for determining the probability dis- 
tribution of each parameter, given the data. PPPM is a 
specific Bayesian implementation. In effect, it is a fancy 
statistical bootstrap with a method for explicitly incor- 
porating parameter fitting errors. The entire approach is 
completely straighforward to implement in a Monte-Carlo 
by brute force, althouth clever and efficient non-brute force 
methods for implementation of the Monte-Carlo exist (van 
Dyk et al. 2001). 

4.1.2. Details of the line search procedure 

We simulated our spectra with parameters determined 
by Bayes formula. These spectra are folded through the in- 
strumental response using XSPEC. The actual line search 
was done by performing a blind search for absorption lines 
between 0.3 and 4.0 keV. The statistic F was calculated 
as, F = (Xc - xD/xl (Bevington (1969); P02) where F is 
obtained from the fit to a continuum model (C) and con- 
tinuum model plus absorption line (U). A blind search for 
spectral features in fake spectra in XSPEC is dependent on 
the initial values and paths (Rutledge et al. 2002). In order 
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to reduce this effect we started the line search from four 
different line energies. We may not find all the significant 
absorption lines when more than four are present in the 
simulated spectrum. However such cases, which entail a 
large F statistic, occur very rarely. To confirm this we per- 
formed Monte-Carlo simulations with an increased number 
of initial starting energies, and found that the number of 
occurrences of F larger than 5 did not change. In search- 
ing for a line from four starting locations we found that 
sometimes XSPEC fit the same feature more than once. 
We corrected for such double-counted events a posteriori. 

The distribution of the statistic F — Ax 2 /xl is shown 
in figure 2. This distribution was determined for each in- 
strument. In formulating this distribution we recall §4.1 
that this corresponds to the distribution of F under the 
null hypothesis, i.e., that there is no line. This is what is 
shown in figure 2, along with the classical F-distribution. 
The value of F for six models is plotted, along with signif- 
icance line. We conclude that there are no grounds for re- 
jecting the null hypothesis (no third line). The results are 
always in the l-3<7 range and only reach the higher level in 
a few isolated instances involving black-body models and 
PN singles data. 

Our results differ from previous work (B03/DL04) for 
two reasons. Firstly our F-statistic is about 50% smaller 
due to our better fits to the continuum. And secondly 
the use of the correct finite distribution for the F-statistic 
drops the significance of an extra absorption line by a huge 
amount, as happened in the likelihood simulations of P02. 
The net result is a more than six order of magnitude re- 
duction in the significance of the strongest line found in 
the PN data. 

Further support for our conclusions is presented below, 
where essentially the same result is arrived at by indepe- 
dent statistical tests. 

4.2. Goodness-of-Fit test for lines in smooth spectra 

This type of test is described in numerous papers, but we 
mainly follow the arguments and notation of Eadie et al. 
(1983). We test the null hypothesis (Ho) that there is no 
deviation of counts with respect to the continuum level 
against the alternative (Hi) that there is such a devia- 
tion. We define a test statistic T, and unless T exceeds 
some critical value T c then we accept H . Under H , T has 
some distribution P(T). T c can then be implicitly defined 
by the formula a = P(T)dT where a is defined by the 
4<T tail of a normal, standard distribution. This procedure 
ensures that the null hypothesis will not be rejected unless 
a 4a burden of proof is met. This test is binary - unless 
T exceeds T c we accept H - This result can be further 
quantified by calculating 1 — C = P(x)dx. 1 — C 
can be converted into units of a using the same method 
as employed for a and provides a rough measure of the 
significance of the absorption feature. 

We need only define our test statistic T and its probabil- 
ity distribution under Hq. We are searching for a deficit of 
counts in a region of the spectrum beyond that predicted 
by our continuum model. We will assume that the region 
of interest is defined by K energy bins around 2.1 keV. 
Define = observed counts in bin i; Sj = rn — b\, bi = 
true continuum in bin i and bi = estimated continuum in 
bin i. We can then formally construct the null hypothesis 



Ho, there is no line, (s, = for all K) against all alter- 
natives defined by H\. This is the classic goodness-of-fit 

test. The natural test statistic is T — Y^li ( ni ~ bi \ . 

t-~ii—l var(rii — bi) 

This can be expressed, under suitable conditions (ap- 
pendix) as 

i=1 b t + o l 

This test statistic is composed of random variables Zi = 

(rii — bi)/ \J (rii — bi) or cquivalently T = zf. We restrict 
this sum to the K bins of interest. Under the null hypoth- 
esis Ho one can show (appendix) that for m not too small 
that the Zi are normal (Gaussian) and standard (mean 
and variance 1). Thus under Ho, T is by construction x 2 
distributed with K dof. If K ^ 30 one can additionally 
show (appendix) that a new variable can be constructed 
y = (T — K)j\J1K and y is asymptotically standard, nor- 
mal. We select a T c (y c ), and if T < T c (y < y c ) we accept 
Ho (there are no counts in deficit of the continuum model) . 
T c is related to a as discussed above. y c is similarly de- 
fined only with the standard, normal distribution replacing 
P(T). 

4.2.1. Results applied to 1E1207 data 

For all the models we considered the observed T was 
much less than T c so we accept Hq - i.e., there is no line. 
Table 3 lists the observed T and T c for various models. 
Note T c changes for each model because of its dependence 
on the continuum level and fit variances within that model. 
Table 3 also indicates 1 — C, to give some feel for the " dis- 
tance" of T from T c in probability space. Because K is 
sufficiently large in this analysis, we formed the random 
variable y from the observed T, as mentioned above. This 
is shown in table 4 and provides a direct comparison of 
how close (far) the measured y is from y c and the " center" 
of the distribution P(y) in units of a. 

The formulation above is correct in searching a region 
where there is a priori reason to expect a line. There is 
a complication in how exactly to determine the K bins of 
interest (other than that they are around 2.1 keV). If the 
spectral feature is narrow and the region of the K bins 
broad, then the feature can be washed out. This is also 
true in the converse case. To be conservative we checked 
varying sized regions of interest around 2.1 keV. The re- 
sults were not significantly changed for regions compara- 
ble to a few times the energy resolution (from less than to 
slightly greater than the widths reported in B03/DL04). 

In such a blind search one must correct for the fact that 
there are N possible regions of interest (N is the energy 
band divided by the width of a single region of interest) in 
which to find a statistical fluctuation of a given level. Very 
roughly this increases 1 — C by the same factor, which is 
<~ 6. A more detailed treatment can determine the ex- 
act factor (Eadie et al. 1983). We have neither corrected 
for this factor in this section nor reduced the significances 
in the next section accordingly. Thus, it should be noted 
that the actual significances are approximately 6 times less 
than indicated in table 3. 

4.3. Significance of the absorption lines by equivalent 

width 
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Fig. 2. — P (F) distribution for PN-single (left) and MOS1 (right) instrument. The dashed line shows the classical F-distribution. 



The final and perhaps most straightforward method for 
deducing an absorption feature's existence is to simply as- 
sume it exists and estimate the significance of the detec- 
tion. In this section we evaluate the significance of the 
residuals around 2.1 keV by fitting a Gaussian absorption 
profile, and then calculate the equivalent width (EW) and 
error bars. This particular absorption line profile gives ex- 
cellent fits to the data (xt ~ 1), so no improvement would 
be had by alternate line shapes. We confirmed other line 
shapes did not statistically significantly affect either the 
EW or the line position. Of critical importance here is 
a proper calculation of the error on the line fit. This is 
a standard non-linear x 2 fitting problem of the type dis- 
cussed by Lampton ct al. (1976) (LMB). We follow their 
approach exactly. The spectrum is fit allowing all the con- 
tinuum parameters to vary freely, and the line depth and 
width are stepped through. The minimum value of y 2 is 
jointly estimated for the line depth and width. A con- 
tour plot of the 68% confidence interval for line depth and 
width is then calculated using the prescription of LMB. 
Using the 68% (t, w) contour we integrate to find the EW 
extrema, and the best fit (r, w) provide the best estima- 
tor of EW. We held the line centroid at its best-fit value 
during this procedure to save computational time. A se- 
ries of test runs indicated that the EW depended only 
weakly on line centroid within the 68% confidence inter- 
val. The contour plots (68% confidence) are shown for the 
PN single and MOS1/MOS2 data in figure 3, for model 
III(BB+BB)*L3. In the case of the MOS data the errors 
in depth and width arc highly correlated, highlighting the 
importance of a joint estimation of the errors using the 
procedure of LMB. 

The results of this analysis are shown in table 4 for all 
the relevant models. In no case is the statistical signifi- 
cance of the line detection larger than <~ 2a. As noted 
in §4.2, this is a firm upper limit, uncorrected for a blind 
search. The strongest residuals are for PN data, and the 
EW for these cases are shown in figure 4. Except for 
the black-body models, the MOS1/2 residuals were essen- 
tially non-existent, so we gave up attempting to fit them. 
In some cases in table 4 the absorption line resided on 
a low temperature black-body component which was so 



weak that no EW could be determined. We note for thor- 
oughness (although we think it obvious) that one must 
calculate the EW using both the low and the high energy 
continuum components, even though the physical model 
may indicate the line only resides on one continuum com- 
ponent. Calculating the EW in this latter case may be 
physically important in interpreting the physics, but in 
assessing the existence of a line only the former procedure 
has statistical significance. 

4.4. Concluding remarks on line search 

We have used three different methods to demonstrate 
there is no spectral line in the 2.1 keV region. A summary 
is presented in table 4. The level of significance of the 
third spectral feature and measures of its significance are 
consistent with each other. In only a few isolated cases 
of black-body continuua does the line rise even to the 3a 
level (and only for the F-test), much less something firmly 
indicative of a line. Since the fourth spectral feature is 
much, much weaker than the third feature, we were un- 
able to analyze it in detail. In §6 below we will present 
plausibility arguments that the statistically insignificant 
spectral residuals in the PN data are due to instrumental 
effects. 

5. PHASE-RESOLVED SPECTRAL ANALYSIS 

While no evidence of a third spectral feature was found 
in the time-averaged data, for thoroughness we decided 
to repeat our analyses for phase-resolved data. First we 
searched for the spin-period using the PN data after cor- 
recting photon arrival times to the barycenter. PN was 
operated in small window mode with time resolution <~ 6 
ms. The MOS data does not have sufficient time resolution 
for spin-period determination. We used XRONOS 5.20 for 
our timing analysis. Using the Epoch- folding method we 
found significant pulsation at P = 424.13076 ± 0.00002 
msec, consistent with the results of B03 and DL04. After 
tagging photons with the best-fit spin-period, we reduced 
spectra in four spin phases. Spectral bins were binned 
to have at least 40 counts in each bin. Figure 5 shows 
a typical phase-resolved spectrum fitted with continuum 
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Fig. 3. — t-w contour plot for a 3rd line (PN, MOS1 and MOS2). The contours were calculated at 68% confidence level. We fit model 
III(BB+BB)*L3 to the data. The crosses show the best-fit parameters at 90% confidence of DL04 (Solid = PN, dashed = MOS1, dotted = 
MOS2). 
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Fig. 4. — Equivalent width [eV] of a 3rd line at 2.1 keV and 68% confidence levels. 



model I(BB+BB) and residuals. We fixed Nh to the fit- 
ted value (1.32 x 10 20 cm 2 from the phase-averaged PN 
spectral analysis). Our results did not change when we let 
Nh vary in spectral fitting. 

Around 2.1 keV, only a few spectral bins deviate from 
the continuum model, and residuals are comparable to or 
smaller than the energy resolution, indicating that they 
are due to statistical fluctuation. This is borne out by 
statistical analysis of the residuals using the approach of 
§4.1. Figure 5 summarizes the results for the statistical 
tests on the phase-resolved data. There are no spectral 



features around 2.1 keV. 

6. PHYSICAL SIGNIFICANCE OF THE RESIDUALS IN THE 
CONTINUUM SPECTRUM 

We have demonstrated that there is no third line in the 
1E1207 spectrum. Thus any attempt to explain the origin 
of < 1-3 a residuals should be considered suspect. We 
have several ground rules in this section. We explicitly 
avoid detailed fitting of line shapes and in-depth searches 
for specific instrumental origins of the residual. It would 
certainly be futile given the strength of the residuals, and 
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Fig. 5. — One of the four phase-resolved spectra and residuals fit by model I (BB+BB) (left) and P(F) distribution for phase-resolved 
spectra (right). The dashed line shows the classical F-distribution. 



we are unlikely to offer deep insights into possibly sub- 
tle instrumental effects. Finally we concentrate on moti- 
vating our suspicions using PN data, since the residuals 
are marginally stronger and thus likely to be more fruitful 
hunting ground. 

6.1. Analysis of residual line centroids 

The analysis of §4.3 yields the centroid of the weak resid- 
uals for each of the models described in §3. A plausible hy- 
pothesis is that the residuals are associated with the Au-M 
edge, which occurs in the PN and MOS instruments. To 
pursue this hypothesis we need to focus particular atten- 
tion on the line centroids associated with model III. If the 
residuals arise in the Au-M5 edge of the instrument, then 
model III clearly corresponds to the only physically cor- 
rect model for the data. So if this hypothesis is true, we 
expect to see residuals around the Au-M5 energy. 

Figure 6 shows the line centroids from fits to the third 
line using model III for PN single data, PN double data 
and combined PN single and double data. We also show 
the B03/DL04 90% confidence interval band for their line 
centroid determination. The dotted line is the theoretical 
position of the Au-M5 edge. The striking feature of this 
plot is that the PN single data, which is without question 
the most reliable data since it rejects split events, occurs 
at the position of the Au-M5 edge. Moreover the poten- 
tial origin of the lower residual energy found in previous 
work is suggested by the fact that the PN doubles (split 
events) residuals occur at lower energy. Our combined PN 
singles/doubles data set is in agreement with the results 
of B03/DL04. Thus we conclude the PN split events have 
reduced the apparent energy of the residuals, while the 
pristine PN singles events are at the expected energy for 
an instrumental effect. One would be tempted to say this 
lower energy in the split events is expected. This would 
get into a discussion about how well split events are re- 
constructed in EPIC data and is far too detailed for us. 
We simply note that the highest quality event data shows 
residuals at the Au-M5 edge and the lowest quality data 
shows it at much lower energy. We also show in figure 6 
the position of the residuals in 3C273, using a fitting pro- 



cedure described in the next section. The 3C273 residuals, 
which are most certainly associated with the Au-M5 edge, 
define the "fit" position of the Au-M5 edge and is a more 
relevant metric than the theoretical position of the edge. 
Nevertheless, the position of the Au-M feature in 3C273 
is consistent at the 90% level with the theoretical Au-M 
energy and it is also consistent with the position of the 
residuals in the PN singles data. A better line shape ap- 
proximation would no doubt move the 3C273 result down- 
wards, but this would require detailed understanding of 
the instrument response. Our results are meant only to 
be suggestive since, after all, the residuals are statistically 
insignificant. 

6.2. Consistency test between PN and MOS data 

In DL04 a large inconsistency (~ ha) between the sig- 
nificance of the lines detected in PN and MOS data was 
reported, despite fit parameters being in statistical agree- 
ment. They explained this as being due to low energy cal- 
ibration uncertainties in the MOS. This observation moti- 
vated us to perform a consistency check on our data set. 
Because we utilized a MOS energy cut slightly higher in 
energy, one would reduce low energy calibration uncer- 
tainties, and thus better consistency between MOS and 
PN data sets. This is what we have found, as described 
below. 

We set up a Monte-Carlo simulation to test whether the 
3rd line residual detected by PN could have appeared in 
MOS data as a deficit with a significance comparable to 
that observed. We used the III(BB+BB)*L3 model fit to 
the data because it showed the largest significance of a 
3rd feature, therefore the largest discrepancy between PN 
and MOS. Hereafter we describe our formalism using PN- 
single and MOS1 case as an example: (1) we fit the model 
III(BB+BB)*L3 to the PN-single data and calculate best- 
fit parameters (2) we fold the model with variance in both 
the continuum and line parameters through MOS1 detec- 
tor response and simulate a spectrum with the same ex- 
posure time and spectral bins as our analysis of the real 
data. (3) we evaluate A% 2 by fitting a simulated spectrum 
with and without a 3rd line. We repeat this procedure 
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Fig. 6. — Line energy of a fitted 3rd line. The shaded region indicates the 90% confidence level found by DL04. All error bars are 68% 
confidence. The dotted line is at Au-M5 edge position. 



1000 times and compute the distribution of A% 2 . Figure 
7 shows the results from the consistency test between PN 
and MOS1 (left) and PN and MOS2 (right). 

We calculated the fraction of simulated A% 2 points 
smaller than or equal to the measured A% 2 (S), i.e. 
/(Ax 2 < S). The results are f(A X 2 < S) = 2.9 x 10~ 2 
(MOS1) and f(A x 2 < S) = 3.1 x 10~ 2 (MOS2). The PN 
and MOS1 results are just consistent at the 1.90<7 level 
and at the 1.86<r level between PN and MOS2. We con- 
clude that the behavior of PN and MOS instruments is 
internally consistent with regards fitting of the lines, as 
we showed previously was the case for the continuum fits. 

6.3. Strength of the residuals in the Au-M region 

There are no statistically significant additional spectral 
features in the 1E1207 spectrum. For thoroughness we 
indicate just how weak the residuals are that we have un- 
covered at the Au-M edge. We do this by considering the 
residuals relative to 3C273, where the same instrument 
residuals are well-established (Kirsch 2004). 

We extracted PATTERN 0-4 events from a 40" radius 
region centered on the source. We used a double power law 
continuum model to fit the PN singles data in the 0.3-6.0 
keV energy range. We did not fit the Au-M region but in- 
stead fit our model III with blackbody spectra for 1E1207 
and double power-law continuum for 3C273, and plot the 
ratio of the residuals to the model spectrum (figure 8). We 
note that the residuals are ~ 6-7% in the quasar data and 
~ 10-15% in the neutron star data. The neutron star data 
has larger error bars. 

We directly calculated the EW of the 2.2 keV residuals 
without a line fit, and find them to be ~ 40 eV. This is 
consistent with the EW we determined by line fitting. A 
similar calculation (as a sanity check) gave us 8 eV for 
the EW residuals at 2.2 keV in 3C273, in agreement with 
measurements of the XMM-Newton team (Kirsch 2004). 
Of course the much higher noise level in the 1E1207 data 
renders the residuals insignificant. In figure 8 we show 
the ratio between the continuum model and the data for 
1E1207. The same plot is shown for 3C273. We see there 
is nothing extraordinary in the NS data. The fluctua- 



tions are uniform through the entire spectrum. In 3C273 
the Au-M edge is readily apparent because of the better 
counting statistics. 

We note that a common error in estimating residual 
strength is to fit absorption lines and a continuum, and 
then to use the resultant continuum without lines as a basis 
for comparison with the residuals. This procedure grossly 
overestimates the continuum, and thus the strength of any 
absorption features with respect to it. The correct proce- 
dure of comparing the dips with respect to the best fit con- 
tinuum markedly reduces the significance of the residuals. 
With this proper procedure our results here are consistent 
with those we obtained in §4. 

7. CONCLUSIONS 

• Utilizing a wider range of physically plausible INS 
models and careful background subtraction we have 
obtained continuum fits for 1E1207 which effec- 
tively rule out the need for additional harmonics of 
the first two spectral absorption features. Our de- 
rived continuum parameters for both PN and MOS 
data are consistent with those found by previous 
investigators. 

• To explore the significance of the weak residuals 
found at higher energies in a few cases, we have im- 
plemented an improved F-statistic analysis for ad- 
ditional weak lines in 1E1207. It corrects for sys- 
tematic overestimation of absorption line strengths 
using the standard F-test, as previously reported in 
the literature. Calculating the correct, finite statis- 
tics reference distribution for the F-statistic we find 
lines at harmonics of the two firmly established lines 
are statistically insignificant. 

• We performed three completely independent statis- 
tical analyses to test for the presence of third and 
fourth spectral features. All tests give consistent 
results and all tests indicate there is no third or 
fourth line. 




Fig. 7. — Results from the consistency test for PN-single data. The grey lines correspond to the measured A% 2 values in MOS1 and MOS2 
data. 
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Fig. 8.— Data/model ratio plot for 1E1207 (black) and 3c273 (red). 



• While there are no statistically significant spectral 
features at higher energies, we explicitly fit the spec- 
tral residuals around 2.1 keV. Using the continuum 
model favored by previous investigators, and for 
which we get excellent continuum fits, we find the 
PN single spectrum residuals occur at the position 
of the instrumental Au-M edge. The residuals ap- 
pear at the same energy as those obtained by fitting 
the known Au-M residuals in 3C273 with the iden- 
tical procedure. We further show that PN doubles 
events yield a much lower energy for the residuals, 
as does the combined PN singles/doubles data set. 
This is highly suggestive of an instrumental origin 
for the insignificant residuals. Indeed in all eight 
physically plausible INS models we considered, the 
residuals always appear at the Au-M edge. 



• We have demonstrated that the strength of the 
weak instrumental residuals in the 1E1207 spec- 
trum are consistent with the stronger instrumental 
residuals observed in the 3C273 spectrum. 
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APPENDIX 

TESTING FOR A SPECTRAL LINE IN A CONTINUUM 

Some important results used in §4.2 are summarized here. The arguments of Eadie et al. (1983) are closely followed. 

Consider searching a region of interest (ROI) consisting of K bins for a deficit of counts over the expected continuum. 
The total spectrum has N bins. Define 6j as the true continuum, b\ as the estimated continuum, and rij as the observed 
counts in bin i. We can define the difference of counts in bin i as s, = n; L — 6j. 

The search for residual in the K bins of the ROI can be formulated as a goodness-of-fit test, and under the null 
hypothesis H there are no residuals, and E(si) = E{rii) — bi = bi — bi ~ bi — 6, = 0, where E is the expectation (mean) 

operator. A reasonable test statistic is T = Y^fLi , where V() is the variance operator. This statistic is a function 

1—1 Virii—bi) x/ 

of the random variables Zi — , ni ~ bi = . 

VV(ni-6i) 

The Zi have the following properties if H is true. E(zi) — 0, and V(rii — bi) = V(nj) + V(6j) — cov(rii,bi) = 
V(bi) + V{bi) — bi + a 2 . covQ is the covariance operator. The trick is that the bi should be estimated using the N — K 
bins outside the ROI, for then rii and b\ arc uncorrelated and their covariance vanishes. The continuum is simply fit 
without utilizing the K bins in the ROI. Also note that V(&») = a 2 can be directly obtained from covariance matrix of 
the continuum fit. Thus under Ho, 

K ( h\2 K 

r = E^^- = 5>< (AD 

i=l b i + °i i=l 

and E(zi) = 0, V(zi) = 1. Thus the £j are standard, normal variables (for Sj not too small), and P(T) is a % 2 distribution 
with if degrees of freedom. 

For K > 30 we can simplify this result by defining a new statistic y — ^= ■ Since the mean and variance of the \ 2 

distribution are K and 2K respectively, y has a standard, normal distribution. Then as described in the text we can test 
Ho by simply calculating y and if y < y c — 4 we accept Ho. 
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n qco+0.074 


n ^ir+ 017 

U.01U_o 040 


U.dbU_ o 6 4 


i?^ [km] a 


^ ri^+o 13 


1 nn+0 19 


-1 r.4-0 19 


~ ~ +ri 40 


1 r>-i +0 54 


n ^rr+0 72 

0-6 7 I .23 


Ex [keV] 


741 + nnnf 




725+"'°"? 

u - ' zu -0.005 


730+^nnr 




728+^nnr 

u. / z,o_ 005 


ri 


Q2+ 05 

u - yz -0.07 


QQ+0.26 

u - yy -o.io 


60+ 06 

u - DU -0.03 


u - oo -0.02 


u - oo -0.01 


u - oo -0.03 


FWHM [keV] 


n ono+0.024 
U.OUO_q 02 g 


297+ 015 
u - zt " -0.022 


9 90+ 030 


n 979+O.O23 
U.z f ^_ .018 


262+ 017 

U.zuz_ 016 


2 64+ - 020 

u - zo ^-0.019 


£2 [keV] 


1 qoo+0.009 
i.oao_ 010 


1 009+0. 008 
i - ooz -0.009 


1 qoo+0.008 
i.ooo_o.008 


1 qni +0.008 
1-O»l_o.008 


1 Q7Q+0.008 
J-.Oi»_o.o08 


1 qoo+0.008 
i.ooy_ 009 




2.9_!. 2 b 


o-58±g:ig 


0.39t°;°2 


u - oo -0.07 


q 7+13 
°- ' -2.0 


u - oy -0.04 


M [keV] 


u - i0 -o.io 


o.i83±8:8?l 


o-i87±8:gg 


iqn+ 027 

U.iyu_ .024 


1 40+ 029 

U.11U_ 052 


n 1 C7+0-028 
U.18/_o.o23 


X 2 /dof 


0.98 


0.98 


0.99 


0.98 


0.98 


0.98 


dof 


422 


422 


422 


422 


422 


422 



Note. — For hydrogen atmosphere models, we fixed neutron star mass and radius to 1.4 M & and 10 
km. kT°° and R°° are insensitive to variation of M and R (Zavlin et al. 1998). 



Note. — All parameters include a roundoff digit except when it would have required displaying a fourth 
digit after the decimal point. 

a For apparent radius, we assumed distance is 2.1 kpc and the error bars do not include uncertainty in 
distance. The distance to the supernova remnant PKS 1209-51 was measured to be <~ 2.1 kpc (Giacani 
ct al. 2000). 

b No upper bound was obtained (see §3.4) 



Table 2 

Best-fit continuum parameters to MOS data 



Instrument MOS1 


MOS1 


MOS2 


MOS2 


Model III (BB+BB) 


III(HA+HA) 


III(BB+BB) 


III(HA+HA) 



N H [10 



kT, 



s 



kT\ 



H 



22 cm 2 ] 
[keV] 
[km] 
[keV] 
i?f [km] 
Et [keV] 

T\ 

FWHM [keV] 
E 2 [keV] 

7-2 

FWHM [keV] 
X 2 /dof 
dof 



n5+ 010 

U - UO -0.04 
1Q7+0.026 
u ' iy ' -0.027 
o 4+1-5 
J -*-1.2 

0.3741°™ 

0.391^7 

7 00+ 033 

u - ' uu -0.031 

90+ - 52 

u - yu -0.22 

u -*°-0.08 

1 4 1O+0-009 
1.418_ 0010 

46+ - 04 

u - q:O -0.05 

236+ 043 

U.ZOO_ 038 

0.85 
139 



1 7fi+ - 035 

U- 1 ' °-0.037 

n 1 1 7+0.025 

u ' ii ' -0.016 

20.31?-? 
41 + ' 11 

-,+0.46 
-0.20 



0.49^ 



n 790+0.013 

U.('28_o.016 
62+ 12 

u - DZ -0.05 

hO.043 
-0.031 
+0.009 
0.009 
+0.05 



0.330 
1.417 

0-39-0.05 

n 1 81 +0-035 
U.l»l_o.030 

0.80 
139 



075+ ' 074 

u - u '°-0.053 

2 04+ 021 

U.ZU1_ .030 
o n+1.6 
°- u -0.8 
+0.056 
-0.039 
0.38 



0.373 



0.37_ .i7 

n 71 O+0. 018 
U -' i8 -0.046 



F0.49 
-0.12 
^0.15 
'-0.06 

1 403+ 010 

1.1UO_0.009 

n 47+0-05 

u -*' -0.05 

2 29+ 035 

u.zz»_ 037 
1.20 



0.68] 
0.39^ 



n 109+0.011 

U - i8z -0.050 
n 190+0. 018 
U.1ZO_ 020 

18l 6 9 

0. 418™ 

ho.59 

-0.17 
-,+0.013 
-0.016 

u - oo -0.06 
a oiq+0.056 
U.Ol»_ 037 

1. ^UZ_ .009 



0.43] 
0.729^ 



+0.04 
-0.04 



137 



0.44 

21 ^+ 034 
u.zio_ 031 

1.18 

137 



Note. - All parameters include a roundoff digit except when it would have 
required displaying a fourth digit after the decimal point. 



Table 3 



GOODNESS-OF-FIT PARAMETERS OF 3RD LINE TEST 



Continuum model 


Parameters 


PN single 


PN double 


MOS1 


MOS2 




T 

± c 


112.8 


71.38 


65.26 


65.26 


III(BB+BB) 


T 


85.90 


43.42 


24.79 


41.03 




1-C 


1.7 x icr 2 


5.4 x 10~ 2 


5.3 x 10" 1 


3.1 x 10~ 2 




T 

± c 


112.8 


71.38 


65.26 


65.26 


III(HA+HA) 


T 


80.45 


42.62 


20.47 


34.68 




1-C 


3.3 x 10~ 2 


6.3 x 10~ 2 


7.7 x 10- 1 


1.2 x 10" 1 



Note. — Model I and II show very similar results as model III. 



Table 4 



Summary statistics for 3rd feature significance 



Instrument 


Continuum model 


P(F) analysis a 


Goodness-of-fit test 


^ b EW analysis a >' 




I (BB*L3+BB) 


1.77 


2.13 


* 




I (BB+BB*L3) 


3.29 


2.13 


2.07 




I (HA*L3+HA) 


2.81 


1.18 


1.69 




I (HA+HA*L3) 


2.84 


1.18 


1.96 


PN-single 


II (BB+BB*L3) 


3.27 


2.22 


2.18 




II (HA+HA*L3) 


2.55 


1.69 


1.85 




III (BB+BB)*L3 


3.26 


2.17 


2.14 




III (HA+HA)*L3 


3.12 


1.83 


1.88 




I (BB*L3+BB) 


0.70 


0.33 


* 




I (BB+BB*L3) 


2.09 


0.33 


1.92 




I (HA*L3+HA) 


0.32 


-0.92 


# 




I (HA+HA*L3) 


0.32 


-0.92 


# 


MOS1 


II (BB+BB*L3) 


2.36 


0.28 


1.82 




II (HA+HA*L3) 


1.02 


-0.83 


# 




III (BB+BB)*L3 


2.55 


-0.07 


1.36 




III (HA+HA)*L3 


1.32 


-0.73 


# 




I (BB*Ui+BB) 


0.69 


1.78 






I (BB+BB*L3) 


2.21 


1.78 


1.65 




I (HA*L3+HA) 


0.27 


1.05 


# 




I (HA+HA*L3) 


-1.43 


1.05 


# 


MOS2 


II (BB+BB*L3) 


2.17 


1.97 


1.68 




II (HA+HA*L3) 


1.14 


1.37 


# 




III (BB+BB)*L3 


2.14 


1.87 


1.87 




III (HA+HA)*L3 


1.11 


1.18 


w 



a The units of significance are sigma. 

b These results have not been corrected for blind search. This will reduce the significance 
of the results (see §4) . 

*Low temperature continuum too weak to determine EW (see §3.4 and §4.3) 

#The upper limit is 1-sigma. 



